Negative thermal expansion of MgB 2 in the superconducting state and anomalous 

behavior of the bulk Griineisen function 
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The thermal expansion coefficient a of MgB2 is revealed to change from positive to negative 
on cooling through the superconducting transition temperature T c . The Griineisen function also 
becomes negative at T c followed by a dramatic increase to large positive values at low temperature. 
The results suggest anomalous coupling between superconducting electrons and low-energy phonons. 
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Superconductivity in the binary compound MgB2 near 
39 K is a fascinating development. Over the last few 
years, scientists have argued that a lattice instability lj 
and/or anomalous phonon behavior might be respon- 
sible for the high transition temperature. Specific atten- 
tion has focused on the E^ g phonon, a bond-stretching 
phonon within the plane of the hexagonal crystal struc- 
ture [3- C3- [H| ■ In-plane tensile strain, induced by lat- 
tice mismatch through thin-film growth, increases the 
superconducting transition temperature T c to 41.8 K; 
this enhancement was attributed to a decrease in the Ei g 
phonon frequency [f| . Phonons can be studied with tech- 
niques such as Raman spectroscopy 0, and heat ca- 
pacity 0, 0, @ . Often neglected in the study of phonons 
is thermal expansion, partly because of the exceptional 
resolution needed to resolve the transition at T c . 

Thermal expansion from powder diffraction measure- 
ments 0, 0| of MgB2 have revealed an anomalous vol- 
ume expansion on cooling below T c . However, high- 
resolution thermal expansion measurements (dilatome- 
try) with a relative sensitivity approximately four orders 
of magnitude better than powder diffraction are required 
for meaningful thermodynamic analysis. Such measure- 
ments of polycrystalline MgB2 were reported but 
discrepancies with the diffraction data, such as the tem- 
perature at which the thermal expansion coefficient a 
changes from positive to negative, are apparent. 

In this Letter, high-resolution thermal expansion mea- 
surements of polycrystalline MgB2 are presented. The re- 
sults reveal a change in sign of a at T c , with negative ther- 
mal expansion below T c ; these data agree with diffraction 
investigations 0, [^} , but offer exceedingly greater detail. 
Analysis of the bulk Griineisen function reveals anoma- 
lous behavior due to dominant low-energy phonon modes. 
The change in sign of a at precisely T c suggests a connec- 
tion between these phonon modes and superconductivity. 

MgB2, synthesized with n B as described previously 
|l2| . was pelletized (diameter = 4.6 mm), placed in a 
boron nitride crucible and heated to 800 °C for 30 min 
at 3 GPa using a cubic multi-anvil press. A very thin 



black layer, impurities from surface reaction with boron 
nitride, was removed, leaving behind a brilliant gold- 
colored MgB2 sample with density 2.56 g/cm 3 (100% of 
theoretical density). Heat capacity was measured with 
a thermal-relaxation technique. T c versus pressure was 
determined inductively (0.12 Oe rms field at 1023 Hz) to 
0.63 GPa using a helium pressure medium; a manganin 
sensor at room temperature served as a manometer |l3j . 

Thermal expansion was measured with a novel capac- 
itive dilatometer cell (Fig. 1) constructed from fused 
quartz [J^, which has among the smallest known ther- 
mal expansions |l5l |. Calibration is realized by measur- 
ing copper at 300 K and 79.3 K, where the length 
change of fused quartz, AL/ L 30 qk, is zero During 
measurements, the sample and cell (in thermal equilib- 
rium) are warmed at 0.2 K/min as capacitance and tem- 
perature are recorded at intervals ranging from 0.1 to 
0.2 K. The data are corrected for temperature-dependent 
changes in capacitance of the empty cell, measured in a 
separate experiment, and the differential expansion of the 
cell/sample using data for quartz [J^. These contribu- 
tions amount to 2.7 % and 4.1%, respectively, of the ca- 
pacitance change due to the thermal expansion of MgB2 
over the temperature range 6 K < T < 300 K; the lat- 
ter correction would be 63 times larger if a conventional 
copper dilatometer cell were used to measure MgB 2 , be- 
cause of copper's large thermal expansion [lfij |. The spec- 
imen was measured along two perpendicular axes (4238 
lim. and 3665 fiin long) to search for preferred crystallo- 
graphic orientation effects (none were detected). 

Linear thermal expansion AL/Lqk is shown in Fig. 
2. These raw data have been corrected for the thermal 
expansion of quartz as described above with no other 
processing. The length decreases with temperature, as 
occurs in most materials, but it expands on cooling be- 
low T c = 38.7 K as illustrated in the inset where the 20 
A scale reveals the absolute length change of the 4238 
fiin long specimen. The distinct change in AL/ Lqk near 
T c is highlighted further in the inset of Fig. 3 where 
angstrom- scale resolution and the continuous (second- 
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FIG. 1: Composed of two L-shaped pieces connected via 100 
/im thick springs, all components of the dilatometer cell are 
made of fused quartz. Cr/Au (100 A/1000 A) 1.6 cm 2 capac- 
itor plates are perpendicular to the page's plane (electrical 
contacts shown); a guard ring (at earth) encircles the low 
plate. A Be-Cu strap (not shown) secures the sample block. 
Expansion of the sample changes the capacitor gap. 
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FIG. 2: Linear thermal expansion AL/Lqk of the 4238 /im 
long specimen. The region near T c is shown in the inset; the 
20 A scale indicates the absolute length change. 

order) nature of the phase transition are evident. The 
measurement was repeated 3 times along the 3665 /jm 
length and 3 times along the 4238 /im length with simi- 
lar results to those shown in Figs. 2 and 3. The overall 
linear expansion {L 300K — L 6K )/L 6K = 119.2(4)xl0~ 5 . 

The thermal expansion coefficient a = 
(l/LsK)dAL/dT is obtained by fitting the data in 
Fig. 2 using Chebyshev polynomials and differentiating; 
overlap between neighboring fit regions assures a smooth 
derivative and caution is exercised near the sharp feature 
at T c . This method is typical for taking derivatives when 
tiny changes between neighboring data points occurs 
|17|. The resulting a values are plotted between and 
50 K in Fig. 3; all six measurement showed similar 
behavior. Note that cv(50 K) is 16.7 times smaller than 
that of copper 0], a consequence of the large Debye 
temperature of MgB2 0. Negative thermal expansion 
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FIG. 3: Inset shows AL/L^k near T c illustrating the change 
in slope, excellent resolution of the data, and continuous na- 
ture of the phase transition. The solid line in the inset is the 
fit that, after differentiation, yields the thermal expansion co- 
efficient q (main panel). Solid lines indicate the jump in a at 
T c . The dashed line below 6 K denotes hypothetical behavior. 

sets in at precisely T c and further inspection reveals a 
crossover at 18.7 K to a > 0. Furthermore, a broad 
peak is evident in a at 9.7 K, below which it decreases. 
According to the third law of thermodynamics, a should 
approach zero at T = 0; this behavior is illustrated by 
the (hypothetical) dashed line in Fig. 3. 

The jump in a at T c (solid lines Fig. 3), Aa = a^- 
as = 6.3(7)xlCP 8 K~ 1 ,was obtained after averaging the 
jumps observed in six measurements; TV represents the 
normal state and S the superconducting state. This value 
agrees with Aa = 5.8xl0~ 8 K _1 reported previously [llj : 
in that report, a < below 10 K, but remained above 
zero for T > 10 K. These differences may be associated 
with: (1) the high density of our specimen, since poly- 
crystalline MgB 2 prepared by the method in Ref.[ll( is 
highly porous with low cohesiveness and elastic modulus 
[13 , (2) the higher T c (38.7 K vs. 37.8 K), or (3) the 
copper thermal expansion cell used in Ref. which 
requires a significantly larger correction to the data than 
our cell (see above). Behavior in AL/Lqk similar to that 
of Fig. 2 was revealed in diffraction measurements pSlllcj. 
but the present data offer far more detail. 

Heat capacity at constant pressure Cp was measured 
on the same sample. A background of the form Cp = 
A\T + A 3 T 3 + A 5 T 5 {A\, A 3 , and A 5 are constants) 
is subtracted, which was obtained by fitting Cp(9 tesla) 
in the range 25 K < T < 50 K; the form of this back- 
ground is identical to that used previously 0- A clear- 
feature is observed at T c = 38.7 K (lower inset of Fig. 4). 
The solid lines denote the step-like behavior expected 
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FIG. 4: The main panel and upper inset show the bulk 
Griieneisen function 7g versus temperature. The lower in- 
set shows heat capacity Cp after subtraction of Cp(9 tesla); 
the solid line indicates the jump expected for a continuous 
(second-order) phase transition. 

for a continuous (second-order) phase transition. Using 
an entropy-conserving construction, the jump ACp = - 
121(1) mJ/mole K is estimated. 
The Ehrenfest relation 



dT c Aa 

1p =3vTc ac- p 



(1) 



provides the connection between the pressure derivative 
of T c , Aa, and ACp (v is the molar volume). It has 
yielded dT c /dP in a gree ment with experiment for other 
superconductors [Tgtl20| . Using Aa = 6.3(7)xl0" 8 K _1 , 
ACp = -121(1) mJ/mole K, and v = 1.74xl(T 5 m 3 /niole 
0, dT c /dP = -1.05(13) K/GPa is obtained. Variation 
in ACp values (66 mJ/mole K to 133 mJ/mole K) @ 
underscores the importance of measuring ACp and Aa 
on the same sample to reliably estimate dT c /dP. 

To test the calculation via Eq. (1), dT c /dP was mea- 
sured on our sample using a purely hydrostatic method. 
Pressure was changed below 55 K and only one tran- 
sition (P = 0.63 GPa) was measured in solid helium. 
T C (P) was determined from the midpoint of the super- 
conducting transition in the real part of the ac magnetic 
susceptibility. The data (Fig. 5) show that T c decreases 
at dT c /dP = -1.07(0.03) K/GPa, similar to prior reports 
[l3j |. This value is in excellent agreement with the value 
predicted by the Ehrenfest relation. This observation 
together with the continuous behavior of AL / L^k elim- 
inates the possibility of a significant lattice instability in 
MgB 2 for 6 K < T < 300 K. 

While negative thermal expansivities are not that un- 
usual 21], the crossover of a from positive to negative, 
at precisely T c , is rare. Tantalum exhibits a slight neg- 
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FIG. 5: In the inset T c versus pressure is plotted. These data 
were determined from the midpoint of the transition in the 
ac magnetic susceptibility Xac, shown in the main panel. 

ative a below T c — 4.4 K; what distinguishes this 
feature in MgB2 is that a is eight times larger above T c , 
and that it becomes strongly negative below T c . Thus, 
the electronic transition into the superconducting state 
has a significant impact on the thermal expansion at and 
below T c as revealed in Figs. 2-3. 

Thermal expansion is a sensitive probe of phonon be- 
havior which is analyzed through the bulk Griineisen 
function 7g = PBV/Cy calculated using data for 
the bulk modulus B and our measurements of the vol- 
ume thermal expansion coefficient (3 — 3a, V and Cp 
(the heat capacity at constant volume Cy ~ Cp It 
is a dimensionless weighted sum 21] of 7^ values (i.e., 7<j 
= J2i c ili/J2i c i) eac h associated with a phonon mode 
of vibrational frequency uji and ji = —dlntOij dlnV . "/g 
has a value of 1.30 near 300 K, rises slightly to 1.35 at 
90 K followed by a gradual decrease. This decrease con- 
tinues (see Fig. 4), until jg crosses zero at T c , reaches a 
minimum of -1.65 at 24.5 K, before increasing to 34.2 at 
5.8 K. Nonmonotonic behavior of the Griineisen function 
is a consequence of strong anharmonic lattice vibrations 
[24| , which have been observed in MgB2 0, 0] . The neg- 
ative as well as the large positive values of "fa rnust be 
associated with anomalous phonon modes that dominate 
the weighted sum due to their respective thermal popu- 
lations; electron-phonon interactions will also play a role. 

In other layered materials, negative a values are gener- 
ally associated with transverse acoustic vibrations, which 
propagate in the layers , leading to negative a-axis 7; 
values. As a result, the change in sign of a and r YG at 
T c , coincident with hardening of the F<2g phonon |4j, is 
likely connected to planar electron-phonon interactions 
associated with the cr-band (in-plane) superconducting 

gap mm. 
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At low temperature, differences in the intralayer bond- 
ing (along the c-axis) may lead to low-energy phonon 
excitations and large c-axis 7, values with high relative 
thermal population that cause the upturn of 7^ below 20 
K; such behavior occurs in other layered materials, albeit 
with significantly smaller values of jg .24]. In MgB2, the 
spectacular increase of 7g coincides with a reduced Cp 
below rs 11 K, attributed to a second superconducting 
gap 0,0- Thus, Cp is diminished because of electronic 
effects due to the second gap, a remains large, and their 
ratio leads to anomalously large 70 values. These com- 
bined considerations suggest that the increase in jg a t 
low temperature may result from out-of-plane electron- 
phonon interactions, associated with the weaker 7r-band 
superconducting gap H^,|2^|, which have a strong impact 
on the c-axis % values. 

The results presented here compliment spectroscopic 
phonon probes, some of which may be surface sensitive, 
by illustrating unusual phonon behavior in two funda- 
mental bulk thermodynamic quantities, the Griineisen 
function and the thermal expansion coefficient. Thermal 
excitations of electrons across the two superconducting 
gaps, and their coupling to low-energy lattice vibrations, 
play a role in the anomalous behavior of 7g and a. 
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